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A generalization of the Drude model is studied. On the one 
hand, the free motion of the particles is allowed to be sub- or 
superdiffusive; on the other hand, the distribution of the time 
delay between collisions is allowed to have a long tail and even 
a non- vanishing first moment. The collision averaged motion 
is either regular diffusive or Levy-flight like. The anomalous 
diffusion coefficients show complex scaling laws. The conduc- 
tivity can be calculated in the diffusive regime. The model is 
of interest for the phenomenological study of electronic trans- 
port in quasicrystals. 

PACS number: 7f.fO.Ca, 05.60.-l-w, 72.fOBg 



Drude's model MM considers a gas of independent 
charged classical particles. Their ballistic motion is per- 
turbed by collisions at random times distributed accord- 
ing to the Poisson process of rare events. At every col- 
lision, a particle exchanges impulsion and energy with a 
scatterer in a random way. The scatterers are supposed 
to be at thermal equilibrium and the temperature of this 
bath partially determines the distribution of the scat- 
tering parameters. The collisions allow the particles to 
dissipate impulsion and energy to the bath, but they also 
force the collision averaged time evolution to be diffusive. 
Placing the system in an exterior electric field, one ob- 
tains the celebrated Drude formula for the conductivity 
in this diffusive regime. 

The first modification of the Drude model considers 
the free motion of particles between collisions. Let us 
suppose that it is characterized by an exponent cr taking 
values in the interval [0, 1] and an anomalous speed v^ 
such that 



x{t) = x{0) + Vat" 



(1) 



where x(t) is the position in euclidien physical space 
at time t and x{0) the initial position at t = 0. Of 
course, such a motion cannot be deduced from a Hamil- 
tonian in the framework of classical mechanics, but it 
does mimic the anomalous quantum diffusion of a parti- 
cle in a disorder or quasiperiodic on-site potential where 
quantum interferences play an important role [g|. This 
has been shown numerically for the Fibonacci and the 
Harper Hamiltonian y as well as for two-dimensional 
models for quasicrystals pi. Furthermore there is an- 
alytical evidence for the Fibonacci Hamiltonian |^. If 
the Fermi level is in a region of localized states, cr — 
and (|l]) describes the hopping from one localized state to 
another. For the Anderson model, one expects diffusive 
quantum motion (cr — 1/2) in three dimensions and at 
low disorder if the Fermi energy is in the band center. 



The particles undergo inelastic collisions at random 
times tn, n € Z, with the moving disorder, notably 
phonons and other particles of the gas. The frozen po- 
tential such as impurity or quasiperiodic potential has 
already been taken into account and leads to the anoma- 
lous free motion (|l|). The direction of the velocity vec- 
tor Va- after collision is supposed to be completely ran- 
dom; the distribution of its modulus shall be given by the 
temperature T of the bath composed by the scatterers. 
Solely the mean CT,a of Iva]"^ will be needed below. In 
the quantum mechanical framework, this constant can be 
calculated explicitly 0. 

Finally, Drude's choice of the exponential law for the 
distribution of the time delay s between collisions is re- 
placed by any probability law, in particular those with 
a long tail. For sake of concreteness, let us consider the 
family of probability measures on [0, cxi) given by 



p^^r{s)ds = 



M 



ds 

(l-|-s/r)i+^ T 



(2) 



with /i > and t > 0. While the exponent /x character- 
izes the tail of the distribution p^^r, t allows to vary the 
(inelastic mean) collision time, namely the first moment 
of P;i,T (whenever it exists, that is /x > 1). Phenomeno- 
logically, increasing the collision time corresponds to low- 
ering the temperature of the bath. Considering the num- 
ber of collisions up to a given time t gives a continuous 
time stochastic process with integer values, a so-called 
counting process. The only counting process which is 
stationary and has independent increments is the Poisson 
process [|j. Hence the counting processes defined by (||) 
do not have this Markov property. Note that, in partic- 
ular, the probability to have no collision after time s + s' 
is bigger than the product of the probabilities to have 
none after time s and none after time s'. Let us finally 
remark that the explicit form of p^^r is of no importance 



for the results below as long as lim sup^ 



,l+M+« 



Pp^,ris) 



is bounded for e < and unbounded for e > 0. Possible 
physical reasons for the choice (||) are not discussed here. 

This accomplishes the presentation and motivation of 
the anomalous Drude model characterized by the param- 
eters fi, T, a and T. The first interesting quantities to 
calculate are the diffusion exponent rj and the (anoma- 
lous) diffusion coefficient Dn{T) defined by means of the 
collision averaged mean square displacement: 



^^,r,a,T{{x{t) ~ x{Q)f) .r^ D^iT)^^ 



(3) 



where 'Ei^T,a,T denotes the mean over all collision times 
and outcomes. Clearly 77 € [0,1]. In order to give a 



precise mathematical meaning to the definition of 77 and 
Dn{T), let us introduce the Laplace transform 



/•oo 

Jo 



fi.T.cr.T 



((f(i)-f(0)f) 



(4) 



Then rj is defined to be half of the infimum over all real 
7 such that Ihns^o S^'^^'Cfi^r, a, riS) = 0. An equivalent 
definition is ^ 

1 ,. Log(E^,.,.,T(x(0-f(0))^) 

77 = - hmsup ^ — - — . (5) 

2 t— 00 Log(t) 

Whenever it exists, the diffusion coefficient is given by 



Dnir) = lim 



i2r? 



lim 5^+2" £^,,,,, J. (,5) 



.5^0 



(6) 



where the equality follows from a Tauberian lemma. The 
Laplace transform C^^r.cr,T can be calculated explicitly: 



^ti.T.a.ri^) 



2aCT,a I^ d,sp^^r{s) Jq dre ^''r 



5rJ2a-l 



1 - ^^ dspf^^r{s)e- 



-Ss 



(7) 



Hence the results can be sumnierized by the phase dia- 
gram given in Figure 1 where 77 and 13,, (t) in the regions 
I to IV are given by 






I '7=2' ^^A'^) = T 

II ,7 = a+i^, Z3,,(t) = r^-1 Ct.. C;,' , 

III ?] = M : 



A,(r)=r2-A'CT^,C-, 



IV 



T] = a , 



A,(r)= Ct^.c;, 



where ( . )^.r denotes the mean with respect to p^^^. and 
C" C'" and C'^ are numerical constants independent of 
T. Note that 77 is well defined and continuous on the criti- 
cal lines separating the four phases, because its definition 
is independent of subdominant logarithmic terms. The 
diffusion coefficient as defined in (o), however, does not 
exist on these lines because of these logarithmic diver- 
gences. Nevertheless, the dependence of D,, on r varies 
continuously across the critical lines. Now follow com- 
ments on each of the four regions. 

The diffusive regime I is remarkable because of the sta- 
bility of diffusion exponent and coefhcient upon change of 
the probability law for the time delays. The dependence 
of Di/2(t) ~ ■y'^"-'^ already appears in pO| , pT| ,p{ . Because 
D^l^if) is linked to the conductivity by an Einstein re- 
lation, this formula is of considerable importance when 







IV 



III 



II 



1 2 1^ 

FIG. 1. Phase diasrain for the anomalous Drude model 



studying electronic transport in quasicrystals. These is- 
sues are further discussed below. In region II, the colli- 
sion rate is not sufficient to force the collision averaged 
motion to be diffusive and hence it stays anomalous su- 
perdiffusive (77 > 1/2). In literature |l^ , p3[ , the motion 
in this regime is generally referred to as a Levy flight and, 
for (T = 1, the above results are known |1^. Levy flights 
can either be obtained from the dynamics of a collisional 
model with long flying times (as above) or alternatively 
with long jumps (as below) (for further references, see 
[|I3|,E4|). The latter approach can be motivated by a max- 
imal entropy principle using Tsallis' generalized entropy 
p4| . Anomalous superdiffusion is of growing importance 
in various flelds of physics such as laser cooling |15[| , ionic 
transport in oxide glasses |lH], non- linear Hamiltonian 
dynamics [[L7[, among others(see the review p3|). In re- 
gion III, the anomalous motion is subdiffusive; the col- 
lision rate imposes the diffusion exponent and only the 
diffusion coefficient depends on the free motion. The case 
CT = is already studied in p8 19|. In region IV, colli- 



sions are so rare that they do not affect the anomalous 
free motion. 

The anomalous Drude model is equivalent to a continu- 
ous time random walk llj,|l^ . This model is defined by a 
probability ipir, s)ds that the walker (or particle) remains 
for a time s at a given site before performing a jump of 
length r between s and s + ds. If ipir, s) = '0(r)'0(s), the 
random walk is said to be separable. For the anomalous 
Drude model, one supposes that the particle stays at a 
given site according to (|l|) and that the jump length is 
determined by (0) according to the length of the waiting 
time. Hence the corresponding continuous time random 
walk is separable only for cr = 0. 

It is interesting to compare the obtained results with a 
corresponding discrete time random walk. In fact, as long 
as /i > 1 (regions I and II), the collision time (s)^^,- is 
finite. It is hence appealing to consider a random walker 
which makes a jump at times N{s)^^r-: A^ G N, with 
spherical symmetric jump probability 



Pt,.T.a,T{.x)d'^x = p^,. 



|1/<T 



-l+l/o- 






drid|a;| , 



(8) 



where |a?| > is the radial variable in spherical coordi- 
nates, dVt the area element on the (rf— l)-sphere and S"^'^ 



the area of the latter. Note that P„ 



(f) 



-i-m/'t 



x{t) = f (0) + Vaf^ + q£ CT.a i 



2a 



(U) 



as |a;| — > oo. The distribution of the position of the walker 
after N jumps is given by the A^-th convolution product 
oi Pp,,T,tT,T- The diffusion exponents ri'{q) and anomalous 
diffusion coefficients D' ,, ^ q > 0, are then introduced 
by 



R<i 



jd p*N 



(x) \x[ 



D',.i,){r)N'i^'^'i^ . (9) 



In region I, /i > 2cr so that the second moment of 
Ptj.,T,a,T exists. Therefore the central limit theorem im- 
plies that the A^-th fold convolution product P*^ „ t^xm) 



with xn = x/yN converges in distribution to a Gaus- 
sian and one verifies that rj' [q) = 1/2 independently of q. 
Of course, this confirms the stability of this regime ad- 
dressed above. Moreover, the ratio of the second moment 
of P^,r.a,T E^nd the coUision time (s)^_t- gives D'-^,r^{T). 
One has D'^,^{t) = Z?i/2(t). 

In region II, Pi_i,^T,a,T{x) ~ |x|~^~'*' with 7 = ^/a such 
that 1 < 7 < 2. Hence Gnedenko's generalized central 
limit theorem Q implies that P*^ ^j^^Xn) with the scal- 
ing xjv = x/N'^ converges in distribution to the spherical 
symmetric Levy law L-y_T,T which is given by its charac- 
teristic function 



Log{L^^r,T{k)) = -C^^T \Tk\' 



(10) 



where Cj^t is a constant depending on 7 and T. Going 
through the details, this allows to determine the diffusion 
exponents given by (9) to be ??'((?) = cr//x for all q < 
n/cr. For q > ij/a the integral in (g) diverges for any 
N. The anomalous diffusion coefficients D' ,, Jt) may 
be calculated as well. 

At first sight, it seems to be paradoxical that ri'{q) 
is defined only for q < fi/a while 77 is defined for q — 
2. In fact, equation (0) gives the stochastic evolution of 
the mean square displacement. It exists for finite times 
because \x{t) — x{0)\ < Ct" if the distribution of the 
Va has a distribution supported in the ball of size C . 
Therefore the diffusion exponent 77 is well defined. On 
the other hand, for the discrete time random walk with 
jump probability P^,r,cr,T, any step can have arbitrarily 
big length so that the moment g = 2 in (0) diverges for 
any A'^. In other words, the distribution of the particle 
positions after a continuous time random walk of time 
N{s)f2.,T is compactly supported while that after N steps 
distributed according to (1q) is slowly decaying. Hence 
there is no paradox at all. 

Note that nevertheless ri'{q) < rj and their values co- 
incide only on the lines separating regime II from the 
others. This is due to the non-Markov property of the 
counting process addressed above which favors long fly- 
ing times and hence big mean square displacement in (g) , 
but has no effect in (m. 

Let us now calculate the conductivity in the diffusive 
regime and verify the Einstein relation. The accelerated 
particle motion in the exterior electric field £ is 



where q is the charge of particle and ct^o- is a constant 
depending on <t and temperature. As M), equation (O) 
does not describe a classical Hamiltonian motion, but 
should be seen as a modelization of the quantum mo- 
tion in an exterior electric field. In fact, if the Hamil- 
tonian quantum dynamics without electric field leads to 
(n[) , then the motion in presence of an electric field satis- 
fies ( [ll| ) in linear approximation in £ 0| . For ballistic free 
motion, (O) gives the usual accelerated motion. If cr = 0, 
(O) means that the hopping between localized states is 
preferably in direction of the electric field. For a = 1/2, 
the motion (O) is ballistic in presence of an electric field 
which allows to calculate the conductivity even if there 
are no collisions. In any other case {a ^ 1/2), collisions 
are necessary as impulsion and energy sink (for a > 1/2) 
or source (for a < 1/2). Strictly speaking, energy dissi- 
pation is only needed when studying thermal transport 
(recall that the electrical current can be calculated in the 
classical Lorentz model although collisions in this model 
are elastic). 

The average speed {v)£„ i^tt (i^ ^^"^ usual sense) of 
a particle evolving between collisions with (O) is given 
by the long time average of ^fj,^T,a,T{x{t) — x{Q))/t or 
equivalently by 

/•C30 

{^Sa^rT = 1™ ^' / di e-**E^,,,,,T(f(t) " f(0)) . 

(12) 

For this limit to exist and to be finite, parameters have 
to be in the diffusive regime I. The calculation gives 

{'^s.cr,t^,r,T ^ 1c-T,crDi/2{T)£/CT,a- If n is the particle 
density, the direct conductivity is thus 



2a-l 2 \^ /m>1 

a = T q ncT.a , , 

■ (s)m,i 



and the Einstein relation reads 



o- = —^ ^i/2(t) • 



Ct,< 



(13) 



(14) 



The anomalous Drude formula (|1^) was obtained in 
|lG,0| and later on in a purely quantum mechanical con- 
text in §,01. 

Using the physical origin of the exponent a as indi- 
cated above, the interpretation of ([l3| ) is the following. 
If cr < 1/2, the particles are trapped by the quantum 
interference phenomena in the on-site potential and, in 
order to have appreciable conductivity, there is a need 
for non-elastic collisions (as in Mott's hopping conduc- 
tivity) . Consequently the conductivity increases with in- 
creasing collision rate (smaller r) and vanishes if there are 
no collisions (r = cxd). On the other hand, if cr > 1/2, 
the conductivity is infinite if there are no collision and 
it decreases with increasing collision rate because the 




FIG. 2. Schematic representation of the temperature be- 
havior of the conductivity as given by the anomalous Drude 
formula 



collisions slow down the free superdifFusive motion. If 
a = 1/2, the conductivity is independent of the collision 
rate. 

The mechanism just described is probably the physical 
origin for the Mooij plot of quasicrystals with differing 
quality. Experiment shows [Eo|, that for those with a 
conductivity bigger than some critical value, the temper- 
ature derivative of the conductivity at low temperature is 
negative, and for those with a conductivity smaller than 
the critical value, the derivative is positive. For the phe- 
nomcnological explanation, one supposes that r ~ j'-a 
for some a > 0. According to (|l^), this leads to the be- 
havior of the conductivity given in Figure 2. Now, the 
higher the quality of the quasicrystal, the lower is the 
diffusion exponent (this expresses the continuous metal- 
insulator transition in quasicrystals). The combination 
of these two facts explains the Mooij plot; the critical 
value corresponds to the materials with a = 1/2. Other 
phenomena such as the inverse Mathiessen rule can be 
qualitatively explained by the anomalous Drude formula 

i- 

In summary, the anomalous Drude model modelizes 
the quantum motion of a particle in an on-site potential 
which further undergoes collisions with phonons or other 
particles by a stochastic process. The diffusion exponent 
and coefficient of the model show complex scaling laws. 
Furthermore the electrical conductivity can be calculated 
in the diffusive regime and the Einstein relation is valid. 
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